To appear in Sept. 10, 2002 ApJ 

Preprint typeset using I£Tf?X style emulateapj v. 14/09/00 



(N 
O 
O 
(N 



OS 



> 
oo 

o 
in 
in 
o 

O 

0,; 

6 
to 

in 
c3 



TIGHTLY CORRELATED X-RAY/Ha EMITTING FILAMENTS IN THE SUPERBUBBLE AND 
LARGE-SCALE SUPERWIND OF NGC 3079 

Gerald Cecil 

Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599-3255 

gerald @ thececils . org 

Joss Bland-Hawthorn 

Anglo-Australian Observatory, Epping, NSW, Australia 
jbh@aao.gov.au 

AND 

SYLVAIN Veilleux 1 
Department of Astronomy, University of Maryland, College Park, MD 20742 
veilleux@astro.umd.edu 
To appear in Sept. 10, 2002 ApJ 

ABSTRACT 

Using Chandra and HST we show that X-ray and Ha filaments that form the 1.3-kpc diameter superbubble 
of NGC 3079 have strikingly similar patterns at ^0."8 resolution. This tight optical line/X-ray match seems to 
arise from cool disk gas that has been driven by the wind, with X-rays being emitted from upstream, stand-off 
bowshocks or by conductive cooling at the cloud/wind interfaces. We find that the soft X-ray plasma has thermal 
and kinetic energies £th ~2 x 1 56 r/° 5 erg and £ke ~5 x 10 54 ?^' 5 erg' where r\ x is the filling factor of the X-ray 
gas and may be small; these are comparable to the energies of the optical line-emitting gas. Hydrodynamical 
simulations reproduce the observations well using large filling factors for both gas phases; assuming otherwise 
leads to serious underestimates of the mass lost in superwinds and therefore their influence within and around 
the host galaxy. X-rays are also seen from the base of the radio counterbubble that is obscured optically by the 
galaxy disk, and from the nucleus (whose spectrum shows the Fe Ka line at 6 keV as well as gas absorbed by a 
moderate neutral hydrogen column). The superbubble is surrounded by a fainter conical halo of X-ray emission 
that fills the area delineated by high angle, Ha-emitting filaments, supporting our previous assertion that these 
filaments form the contact discontinuity /shock between galaxy gas and shocked wind. This X-ray emission is not 
significantly edge brightened, indicating a partially filled volume of warm gas within the shocked wind, not a shell 
of conductively heated gas. About 40" (3 kpc) above the galaxy disk an X-ray arc may partially close beyond the 
bubble, but the north-east quadrant remains open at the surface brightness attained by Chandra, consistent with 
the notion that the superwind has broken out into at least the galaxy halo. 
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1. INTRODUCTION 
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Winds are believed to play a key role in determining the 
equation of state of the intergalactic medium and the long-term 
evolution of galaxies, yet direct proof is elusive. Wind ener- 
getics have proved very difficult to constrain because of un- 
certain gas filling factors in the different phases, but evidence 
is growing that they may be far larger than inferred from op- 
tical spectra of the shock or photoionized nebulae, and radio 
observations of supernovas and their remnants. Winds can in- 
fluence galaxy evolution in many ways (Heckman 2002). They 
are important redistributors of metals by radiation on dust or 
by over-pressured gas from supernova heating. In the Galaxy, 
enrichment of the halo and outer disk through powerful winds 
may account for the observed abundances of the thick disk and 
globular clusters (Freeman & Bland-Hawthorn 2002). More- 
over, winds can regulate the growth of the central bulge (Carl- 
beg 1999), and alter the evolution of satellites of large galaxies, 
with either galaxy hosting the outflow (Irwin et al. 1987, for 
example). Winds may mix astrated gas into intergalactic clouds 
(Tripp et al. 2002), and may clear paths through dusty veils for 
AGN/starbursts to reionize the universe. 



A superwind is driven by the pressure of thermalized super- 
nova or AGN ejecta. As the wind expands and cools adiabat- 
ically, it accelerates to several thousand km s" 1 . Expanding 
along the pressure gradient in the galaxy, the wind can over-run, 
crush, and heat disk and halo clouds to X-ray emitting temper- 
atures, then entrain them into the flow. X-rays may also arise 
from stand-off bow shocks in the wind upstream of the opti- 
cally emitting clouds. Thus, soft X-rays may point to powerful 
shocks. However, they can also be emitted at the conductive 
interface between hot shocked wind and cooler ISM swept up 
by a bubble. Whether gas blows out of the galaxy depends on 
hydrodynamical details and halo drag (Strickland 2001). 

To understand this ubiquitous phenomenon, we are studying 
energetic outflows in nearby galaxies. Here, increased spatial 
resolution of many ISM phases may uncover evidence for pre- 
vious outflows and energetic stellar populations, constraining 
the duty-cycle of ejection and illuminating the feedback role 
of fast and slow winds on supernovae (Efstathiou 2000). One 
of the clearest examples of a nearby superwind-blown structure 
is the prominent bubble (Ford et al. 1986) that protrudes from 
the nucleus of the nearly edge-on SB(s)c galaxy NGC 3079 
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(17.3±1.5Mpc, so l" = 84pc). Rising 1.3 kpc (15") above the 
disk, the bubble is the optically emitting half of a double-lobed 
radio structure (Duric & Seaquist 1988), optical line emission 
from the other lobe being extinguished by the galaxy disk. Out- 
flow kinetic energies of up to 5 x 10 55 erg are implied by the 
optical line emission (Veilleux et al. 1994). Because of consid- 
erable extinction, the nuclear power source is manifest primar- 
ily by reprocessed radiation. How the superbubble is powered 
is thus uncertain, with both AGN (Hawarden et al. 1995) and 
starburst (Sofue & Vogler 2001) sources promoted. 

The total energy required corresponds to the effects of 10 4 /e 
supernovae, where e is the thermalization efficiency (Cheva- 
lier & Clegg 1985). A particular obstacle to achieving realistic 
simulations of the wind is that even a rough estimate of e is 
controversial (Strickland 2001): it can be <§C 10% in dense gas 
(Thornton et al. 1998) where the energy is radiated away, but 
can exceed 30% if the supernovae occur in multi-phase gas. 

In Cecil et al. (2001, CBVF hereafter), we analyzed 
high-resolution images of the superbubble obtained with 
HST/WFPC2, see Fig. la. The Ha+[N II]A6583 emission- 
line image shows many plausibly windblown features, includ- 
ing: 1) linear filaments rising several kpc above the disk to form 
an "X" that we have interpreted to be the contact discontinu- 
ity/ISM shock between wind and galaxy "thick disk'Vhalo gas, 
2) four towers of intertwined filaments that together form the 
optical emitting walls of the superbubble and that all disperse 
at the same height above the disk. The latter pattern, together 
with 3) a break in the velocity field at this point, argue for a vor- 
tex in the top half of the bubble caused by fluid instabilities, as 
predicted from numerical simulations of windblown bubbles. 
4) Line of sight velocities that increase with height above the 
disk to reach 1000 km s -1 relative to systemic at the vortex, de- 
creasing thereafter, 5) a circumnuclear disk region "scoured" of 
dust and gas whose boundary joins the linear filaments 1). 

X-rays could arise when the Ha filaments moving at the 
speeds seen in the superbubble collide with stationary gas, and 
indeed Pietsch, Trinchieri, & Vogler (1998, PTV hereafter) 
found that NGC 3079 is ten times more X-ray luminous than 
galaxies of similar optical luminosity. They used the ROSAT 
HRI to image X-rays across the superbubble with integral 0.1- 
2.4 keV luminosity 10 40 erg s" 1 , and also used the ROSAT PSPC 
to detect the softest X-rays out to 2.'5 east of the nucleus and 
6' west. No spectrum could be isolated from the superbub- 
ble because the HRI was energy insensitive while the PSPC 
lacked spatial resolution to connect optical and X-ray emitting 
filaments. Now a Chandra X-ray image is available, and the 
striking correspondence it reveals and the implications on soft 
X-ray emission from superwind bubbles are the subjects of this 
paper. In §2 we introduce the X-ray image, and analyze it in 
§3. In §4 we discuss the dynamics, and summarize in §5. 

2. OBSERVATIONS AND REDUCTIONS 

We analyzed the archived 26.9 ksec exposure (P.I. D. Strick- 
land) obtained with the ACIS-S3 detector on the Chandra X-ray 
Observatory; the other ACIS-S chips miss the galaxy so are not 
discussed in this Letter. We extracted a 2' square region cen- 
tered on the superbubble (the aim point of the observation); well 
exposed point sources are round in this region, confirming the 
excellent spacecraft aspect solution. A few background "flares" 
were removed by using the ciao2.2.1 lightcurve task, reducing 
the effective exposure to 26.0 ksec. We formed several images 
with different smoothing, using constant and adaptive kernel 



sizes. For study of the nucleus and superbubble we obtained 
best results using a fixed kernel of 0."8 FWHM. Elsewhere, we 
mapped fainter emission using a kernel of 2" FWHM. 

Fig. lb registers the X-ray image with the HST line image 
discussed in CBVF. As others have found in recent studies, the 
"world coordinates" derived from calibrations of both space- 
craft agree within 0"7, which is adequate given the smoothing 
required for this fairly shallow Chandra exposure. We con- 
firmed registration to this accuracy using centroids of a few 
compact galaxies visible in both wavebands. 

To form spectra (Fig. 2), we subtracted the background 
summed from regions off the galaxy to successfully remove 
spectral artifacts, and mapped pulse heights to energies us- 
ing the 32 x 32 response values at the bubble position in 
CALDB file acisD2000-01-29fef_phaN0001. After subtract- 
ing the background, 1757 X-ray events were detected within a 
27x22" ellipse centered on the superbubble (Fig. 2b) and ex- 
cluding a 5. "9 x 3 "9 ellipse that spans the nucleus and jet and 
that emitted 309 events (Fig. 2a). We used the Sherpa pack- 
age within the ciao 2.2.1 software to quantify the spectra and to 
integrate fluxes. We fit the spectrum of the bubble with single 
temperature, collisional-ionization-equilibrium models such as 
mekal and external photoelectric absorption, allowing neutral 
absorption, gas temperature, and interpolated gas abundances to 
vary over the fitting energy interval 0.15- 2.5 keV; abundances 
relative to He were fixed to the Anders & Grevesse (1989) ra- 
tios. For the nucleus+jet spectrum, we also included a single 
Gaussian at the Fe emission. 

Fit results, including derived X-ray fluxes and luminosities, 
are summarized in Table 1. We find that iron and a-element 
abundances are <~ 15 ±5% solar across the bubble (<~ 8 times 
larger than Strickland et al. (2002) derived for the diffuse halo 
of NGC 253), and about half solar near the nucleus. The lack 
of independent constraints on element abundances from opti- 
cal spectra added uncertainty on the fitted emission "integral" 
EI = n e nHrj x V, where r/ x is the X-ray volume filling factor. It 
is possible that depletion onto dust explains the low iron abun- 
dances in both regions, but this would not explain depleted a- 
elements if confirmed by deeper spectra. Calibration uncertain- 
ties (see §2.1.4 in Strickland et al. (2002), notably the effective 
area of chip S3 is uncertain below 0.5 ke V) and low counts from 
both regions limited possibilities for detailed spectral analysis; 
in particular, the necessary spectral binning evident in Fig. 2 
prevented us from using the intensity ratios of successive oxy- 
gen ionization states ( VII/VIII) to check the overall temperature 
solution hence the questionably low fitting abundance of the 
bubble. Low counts also prohibited multi-temperature models. 
We did not explore the addition of a non-thermal component 
to increase gas abundances in the thermal gas although this is 
clearly plausible physically not only near the nucleus (where it 
might be the AGN), but also in the superbubble (where it might 
be synchrotron emission from cosmic rays accelerated locally 
in shocks). 

3. PHYSICAL PROPERTIES OF REGIONS EMITTING X-RAYS 

We now discuss the X-ray emission, proceeding from small- 
est to largest scale. We do not discuss the considerable number 
of X-ray sources evident in the galaxy disk beyond 20" radius. 

3.1. Extended Nucleus 

The nucleus is obscured optically but prominent to Chan- 
dra. X-ray flux is elongated along the axis of the radio jet 
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(§3.2). The nominal registration puts the nucleus 0."4 south- 
west of the X-ray flux centroid, of order our registration un- 
certainty. Indeed our spectral model shows that flux near the 
nucleus is absorbed <~ 10 x more than Ptak (1997) estimated 
from the combined ASCA+ROSAT PSPC spectrum; we could 
not fit the Chandra spectrum with his prefered absorption. Fig. 
2a shows relatively strong Fe Ka near 6 keV, providing unam- 
biguous evidence for the AGN that generally excites this line by 
fluorescence of the cold molecular gas evident nearby (Sofue et 
al. 2001). Sparse counts in the adjacent "continuum" prevented 
us from establishing the equivalent width of this line to compare 
to its strength in other AGN. 

There is no evidence in the spectrum for a hard-component 
due either to the AGN or to unrelated hard-spectrum X-ray bi- 
naries down to a count rate of 0.004 per second integrated over 
the interval 3-7 ke V. 

3.2. Optical Extension of the VLBA-Scale Jet 

Panels A of Fig. lc show anticorrelated X-ray and high- 
velocity optical line filaments associated with the VLBI-scale 
jet (CBVF). The main optical filament (top one in panel A) is 
displaced from an X-ray source (top oval in the top panel A) that 
does not show the hard spectrum of an unrelated X-ray binary 
down to the aforementioned count rate. The X-rays are likely 
not inverse Compton scattered because there is no strong IR 
emission or collimated radio flux along this axis at the scale of 
the Ha/X-ray emission. Instead, the thermal spectrum that we 
infer with kT <~ 0.5 keV is consistent with that expected for a 
jet cocoon shock with average expansion or entrainment veloc- 
ity V sl = y/T s x ts /3.2x 10 5 = 600 km s" 1 , where V s7 = V shock / 1 00 
km s" 1 , and x ts = 2.3 is the concentration of all particles relative 
to hydrogen nuclei when the preshock medium is fully ionized 
and has 10% helium abundance by number (Hollenbach & Mc- 
Kee 1979). This velocity is consistent with measured emission- 
line half-widths in the region, see Figs. 1 lb & c of CBVF. 

3.3. Bubble "Towers" 

Fig. 1 shows that, within the superbubble, optical line and 
X-ray filament networks coincide within our registration uncer- 
tainty (< 0."7 or 60 pc). In both wavebands, the two northern- 
most towers (numbers 3 and 4 in Fig. la) are less diffuse than 
the southern pair. All towers extend down to the galaxy nu- 
clear region. CBVF show (see also Fig. 3) that radio emission 
from the eastern superbubble peaks at larger galactic radii than 
the optical line and X-ray emission, arguing against an inverse 
Compton origin for the X-rays. 

Because of the similar appearance of the superbubble in the 
two wavebands, we bounded the volume of warm gas by assum- 
ing that it envelopes the optical-line emitting filaments. The ar- 
eas of the optical filaments are often upper limits, constrained 
by our inability to resolve many of them in the [N II]+Ha line 
HSTAVFPC2 image. From Table 3 of CBVF we summed the 
volume of those that we now see are associated with X-rays, 
then increased this by a conservative factor of 16 in cross- 
sectional area (assuming cylindrical filaments) to reach the res- 
olution of the Chandra observations, yielding total emitting 
volume V < 10 8 pc 3 . We used the luminosities derived from 
our X-ray spectral fits to constrain the average gas densities 
and other physical properties of the brightest parts of the su- 
perbubble; given the small number of X-ray counts, the derived 
physical parameters are at best order of magnitude estimates. 



The luminosity of the X-ray emitting gas 

L(T) = V x n ex n px A(T) erg s" 1 , 

where n ex (n px ) is the electron (proton) density of the hot gas, 
and A(r) is the cooling function that is approximated near 
r 6 =10 6 Kby 

A(T) = 1.6 x 10- 22 r f -°' 7 + 2.3 x 10- 2 V7^erg cm 3 s" 1 

(McCray 1987). Following Strickland et al. (2002), we scale by 
factor 1Z = 0.15/Z true to show the dependence on abundances; 
in the following j3 = \fr) x 1Z. The X-ray electron density is 
n ex = -y/EI lZ/Vrj x = 18 (3/rj x cm" 3 , while the thermal pressure 
is Pjn/k= 1.7 x 10 & (3/r] x K cm -3 . The mass of the warm X-ray 
gas is M x = 4.5 x 10 7 /3 M Q . The time to cool to 10 4 K at so- 
lar abundance is t c w 3.2 x 10 3 /« c r 3/2 s = 0.15/3 Myr (Suther- 
land & Dopita 1993). The mass flow rate M aow ~ M x v/z = 
296/3viooo M Q yr" 1 , where z = 1200 pc is the height of the bub- 
ble, and viooo is the mass weighted, average outflow velocity in 
units of 10 3 km s" 1 (typical of measured values for the optical 
filaments, CBVF). The soft X-ray plasma has thermal and ki- 
netic energies Ej H = 2 x 10 56 /3 erg and E KE = 5 x 10 54 /3v 2 000 erg. 
In comparison, we found (CBVF) that the optical-line emitting 
gas has E TH = 3 x 10 54 /n«, erg and £ke = 2 x 10 54 f^- 5 v 2 o erg 
where rj v and n ev are the filling factor and electron density of 
that gas phase. 

The mass of both the visual line- and X-ray emitting gas 
can be written as X(l cm" 3 /n e ) M , because the emissivity 
in both cases is proportional to n p n e . From the HST line im- 
age, we derived (Table 3 of CBVF) that n ev > 4.3?7~ ' 5 and 
M v = 1 .3 x 10 6 ?7° 5 M Q . If we assume that the two phases are in 
local pressure equilibrium, we find the relative values of elec- 
tron density (n ev /n ex = 0.2), and mass (M v /M x = 0.03). 

3.4. X-rays From the Radio Counter-bubble 

Fig. 3 maps X-rays over a larger scale and compares to the 
radio continuum. A pair of faint ridges of X-rays emerge from 
the nucleus and trend westward, toward the bottom. They have 
the same opening angle as the bottom half of the bubble walls 
(forming a cone), suggesting that they trace the walls of the base 
of the radio counterbubble. If of similar surface brightness to 
the eastern bubble, they are attenuated about threefold in count 
rate by absorption in the galaxy, so are too faint for meaningful 
spectral analysis. 

3.5. X-ray Emission Near the Disk, Outside the Superbubble 

Fig. 3 shows that diffuse X-ray emission surrounds the nu- 
cleus out to 25" radius. Emission is more extensive east of 
the nucleus, the disk nearside, than to the west. X-rays fill a 
cone with opening angle ^90° as height increases above the 
disk, and are more concentrated on the south side of the nu- 
cleus where the cone edge coincides with inclined linear Ha 
filaments. X-ray emission on the north side of the nucleus is 
patchier, but may also be bounded in the cone sketched in Fig. 
3. The cone apex is not at the galaxy nucleus. Instead, the cone 
reaches the disk while still open, then becomes "two horned" 
west of the nucleus as X-rays trace the radio counterbubble 
(discussed previously). Again, X-ray emission is too faint for 
meaningful spectral constraints. 

3.6. X-ray Emission Above the Galaxy Disk 

We have not analyzed the statistical significance of other X- 
ray point sources. However, one of the registration sources in 
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both wavebands is a small galaxy in the HST I-band image 32" 
east of the nucleus (source "A" in Fig. 3). An arc of 4 sources, 
each about half the X-ray brightness of this galaxy but among 
the brighter "background" sources, lies ^20" south and may 
form the eastern extension of the southern edge of the X-ray 
cone. It curves toward the cone axis as height increases from 
30 to 40". 

4. DISCUSSION 
4. 1 . Outer Wind/ISM Shock 

Proceeding inward from the largest scale, the region of dif- 
fuse X-ray emission is bounded by the optically emitting X- 
pattern of filaments, which CBVF showed have the shape ex- 
pected (Schiano 1986) for the contact discontinuity /ISM shock 
associated with lateral stagnation of the superwind in the galaxy 
"thick disk'Vhalo; the X-ray emission is not significantly edge 
brightened, suggesting a partially filled volume of warm gas 
within the shocked wind, not a shell of conductively heated gas. 
CBVF's models appropriate for steady-state wind flow, and re- 
produced as dotted and dashed white lines in Fig. 3, remain 
open at the top. However, Sofue & Vogler (2001) have mod- 
eled an adiabatic "bipolar hypershell" that arises from a point, 
impulsive energy release on a timescale much shorter than the 
several Myr expansion time of the superbubble. These struc- 
tures close at the top. The appearance of their models is sensi- 
tive to the density distributions of disk, halo, and the density of 
the extragalactic medium, but resemble one another for roughly 
the first 10 Myrs. The outer shell expands at 200 - 300 km s" 1 to 
shock-heat gas to 0.2 keV (so would be very "soft" in a "hard- 
ness" map derived from a deeper Chandra exposure), consis- 
tent with the velocities in our Fabry-Perot spectra (Veilleux et 
al. 1994). Of course, identical shock envelopes arise for con- 
stant values of Eq / po, where Eq is the energy of the explosive 
event and po the ambient gas density. Strickland et al. (2002) 
considered this model to interpret their Chandra image of the 
starburst galaxy NGC 253, rejecting it for one that describes 
optical emission from a wind-blown shell (Weaver et al. 1977), 
the model that (Veilleux et al. 1994) have applied to NGC 3079. 

Sofue & Vogler (2001) compared their models to PTV's 
PSPC image of NGC 3079, noting double-horned features on 
both sides of the galaxy at radii of ~400". This is much larger 
than the scale we are discussing, suggesting either that there 
have been several outbursts over the last 100 Myrs or that the 
larger scale double-horns seen by the ROSAT PSPC are actually 
spurious overlaps of unrelated background objects. 

As discussed in §3.3, what powers the superbubble remains 
ambiguous despite detection of Fe line emission (Fig. 2). As 
emphasized by Maloney (1999), it is extremely difficult to dis- 
tinguish the action of a starburst from an AGN based solely on 
the ionizing flux. Deeper images with Chandra and especially 
XMM/Newton could establish if the contact discontinuity/ISM 
shock closes to form symmetric bubbles. Open morphology 
would argue for a steady wind rather than impulsive starburst, 
supporting the interpretation of Irwin et al. (1987) that H I in the 
companion galaxy NGC 3073, 10' west, trails away from NGC 
3079 because of what would have to be an extremely energetic 
wind. A steady wind would favor an AGN power source, be- 
cause repetitive starbursts would make numerous, less massive 
stars that are not evident (cf. those in M82, de Grijs, O'Connell, 
& Gallagher 2001). 

4.2. Filament Correlation Within the Superbubble 



The organization of both filament systems into vertical tow- 
ers attached to the galaxy disk excludes the possibility that the 
optical filaments are halo clouds caught in the wind. For the su- 
perwind in NGC 253, Strickland et al. (2002) dismiss the pos- 
sibility that diffuse X-ray and filamentary Ha emission both 
come from a cooling, volume filling wind; their argument is 
also valid for NGC 3079. 

That the towers emit X-rays along their entire optical length 
is an important constraint, because the radial velocities of the 
optically emitting filaments show a linear gradient from 200 to 
> 1000 km s" 1 with height above the disk (Veilleux et al. 1994) 
until they overturn in a vortex in the top third of the super- 
bubble (CBVF). The ubiquitous X-ray emission suggests that 
the superwind has driven cool disk gas into the halo, with X- 
rays being emitted either as upstream, stand-off bow shocks or 
by cooling at cloud/wind conductive interfaces. In either case, 
the X-ray emitting skin is thin and would plausibly coincide at 
our resolution with the optical line emitting filaments, e.g. for 
preshock density «o cm" 3 , the separation would be 1 x V s 4 500 /ho 
arcsec scaled for shock speed in units of 500 km s" 1 . In this 
picture, the wind strips gas from the walls of the dense disk cav- 
ity, and this gas then fragments into disintegrating streams by 
Kelvin-Helmholtz instabilities. The fragments gradually dis- 
perse, and eventually become Rayleigh-Taylor unstable in the 
top third of the superbubble where they pass through the vortex 
and begin to drip toward the disk. The filaments are embedded 
in the wind that flows past them to reach into the halo. 

The soft X-rays may arise as optical line-emitting clouds 
evaporate into a region of unseen, hot (10 8 K) X-ray gas that 
fills the superbubble. Because the cloud column density is 
surely < 10 24 cm" 2 , we can use the saturated thermal conduc- 
tive rates determined by Krolik, McKee, & Tarter (1981). The 
approximate cloud lifetime for saturated thermal evaporation 
without an inhibiting magnetic field is 

r TE «0.3Myr/V 2 7 3 /6 r 4 1/6 / J P9, 
where Pg is the thermal cloud pressure in 10 9 erg cm" 3 inT = 
10 7 K cm" 3 ), T4 is the cloud temperature in 10 4 K, and N23 
is the cloud column density in 10 23 cm" 2 that is unfortunately 
unconstrained by existing optical spectra of the superbubble. 
However, because clouds reach high into the galaxy halo from 
the disk, we can set a lower limit on their column densities 
by equating tte to the cloud crossing time z/viooo ~ 1.2 Myr, 
obtaining N23 ~ 0.87^' ' P^~' ' . The visible clouds that have 
survived conductive evaporation are therefore those with col- 
umn densities > 10 22 cm" 3 and hence large sizes (if cylindrical, 
r > 0.3 pc) and total masses M > 10 M Q , consistent with the 
ionized masses determined in CBVF. 

4.3. Wind Energetics 

The detection of Fe K line emission addresses the question 
of whether the AGN is directly responsible for the superbubble. 
Extended filamentation is seen in both Seyfert galaxies (where 
circumnuclear starbursts are commonly observed), and in star- 
burst galaxies that do not show strong nuclear X-ray sources. 
Starburst galaxies commonly show an ionization gradient in the 
extended filaments, in the sense that inner filaments appear to 
be photoionized while the outer ones are more shock-like (Mar- 
tin 1997; Shopbell & Bland-Hawthorn 1998, for example), al- 
though Veilleux & Rupke (2002) have found a starburst-driven 
filament system that is shock-like throughout. From morphol- 
ogy alone, it is difficult to distinguish a shock from a photoion- 
ization process: hydrodynamics can shape filaments that are 
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then photoionized by the central starburst. In the case of M82, a 
well defined ionization cone is mapped in e.g. [N II]/Ha (Shop- 
bell & Bland-Hawthorn 1998)). The optical cone appears to be 
mostly photoionized from stars whereas the X-ray cone appears 
to be shock excited. In NGC 1068, where winds and radia- 
tion are both known to be important, the AGN radiation field is 
thought to explain both the optical and X-ray emission, and this 
object appears to have an X-ray photoionization cone (Young, 
Wilson, & Shopbell 2001). On UV energetics alone, it is also 
impossible to distinguish between AGN or starburst (Maloney 
1999, for example). A gas temperature derived from both the 
X-ray and the optical emission is crucial to establishing the en- 
ergy source. The mean energy of an ionizing photon in pho- 
toionized gas is much smaller than that required to ionize the 
medium collisionally, hence shock heated gas is considerably 
hotter. Suitable optical data do not yet exist for the superbubble 
of NGC 3079. 

True wind energetics are elusive because they depend on 
both the supernova source energetics and the thermalization ef- 
ficiency e. Before Chandra, X-ray observations of starburst 
cores had lacked spatial resolution to resolve the contributions 
of X-ray binaries. Now, ACIS on Chandra can search for hot 
gas in the starburst core, and Griffiths et al. (2000) argue that 
they have detected it in M82. Summing the energy in a wind 
like that in M82 with best estimates for the filling factors in var- 
ious gas phases, one finds that it is still substantially below the 
energy derived using the radio-determined supernova rate and 
a moderate to high (20-50%) thermalization efficiency. More- 
over, Chevalier & Fransson (2001) argue that the supernova rate 
in M82 may be substantially higher than inferred from radio 
counts. These suggest that we are missing much of the outflow 
energy, in agreement with multi-phase hydrodynamical models 
(Suchkov et al. 1994; Strickland & Stevens 2000). 

The problem is that the outflowing gas is almost all entrained 
disk and halo gas, not wind material. As evidence for this, 
the outer filaments in M82 (and NGC 1482, Veilleux & Rupke 
2002) share the rotation of the galaxy disk, and the pattern of 
magnetic fields and gas motions in NGC 3079 (CB VF) strongly 
suggest that the top of the superbubble is forming a vortex. Sec- 
ond, hot gas exists in superwinds (Heckman et al. 2001) and 
appears to be moving faster than the cooler gas, inconsistent 
with standard superbubble models (Weaver et al. 1977) but re- 
sembling the pattern seen in hydrodynamical wind simulations 
(Suchkov et al. 1994; Strickland & Stevens 2000). 

Evidence for extensive supernova remnants in the disk of 
NGC 3079 is ambiguous (see §4.3 of CBVF). We thus address 
wind energetics from the constraints derived previously on the 
efficiency of thermalization. The energetics of the soft X-ray 
component in the superbubble that we constrained in §3.3 are 
comparable to those determined by Strickland et al. (2002) for 
the diffuse wind in NGC 253. This arises even though the X-ray 



emitting filaments in NGC 3079 may have much higher density. 

5. CONCLUSIONS 

Chandra reveals a remarkably tight correspondence (<0."8, 
65 pc) between optical line and X-ray emitting filaments in the 
superbubble of NGC 3079. The bubble is embedded in a coni- 
cal distribution of more diffuse X-rays. West of the nucleus, this 
emission extends into the optically obscured, radio counterbub- 
ble. East of the nucleus (i.e. on the near side of the galaxy) a 
funnel is also delineated by extremely narrow optical line emit- 
ting filaments that rise several kpc above the galaxy disk, sup- 
porting our earlier contention that they form the contact discon- 
tinuity/halo shock of the wind. These filaments connect back to 
the disk at one end and at the other, perhaps, to an X-ray ridge 
that seems to curve back toward the cone axis. However, at 
present X-ray sensitivity, the loop does not close completely, so 
the wind may escape the galaxy. A deeper Chandra exposure 
would confirm this X-ray extension, and would double spatial 
resolution of the superbubble to tighten the volume of the emis- 
sion integral hence constraints on derived physical parameters. 

Even without a complete understanding of the physical pro- 
cess, the close association between the Chandra and HST emis- 
sion is important for several reasons. First, powerful winds 
must be involved either directly or indirectly in producing this 
filamentary complex. Second, high resolution joint Chan- 
dra/HST studies of the filaments in this galaxy, M82 (Grif- 
fiths et al. 2000), NGC 253 (Strickland, Heckman, Weaver, & 
Dahlem 2000; Strickland et al. 2002), and NGC 1068 (Young, 
Wilson, & Shopbell 2001) show that both the X-ray and opti- 
cal emission is highly clumped, with resolved structures seen 
down to our resolutions of 65 pc (Chandra) and 15 pc (HST). 
The evidence for entrainment in these winds is compelling (e.g. 
CBVF), and we suspect that most of the superwind mass resides 
in the filamentary complex. We agree with the assessment of 
Strickland (2001) that the smallest cloud scales in simulations 
must be very well sampled to treat cloud entrainment and ac- 
celeration correctly (Klein, McKee, & Colella 1994). The reso- 
lution of present wind models (Suchkov et al. 1994; Strickland 
& Stevens 2000) is ten times too large to match the observed 
clumping scales, arguing for a new generation of hydrodynam- 
ical models to establish testable predictions on cloud internal 
properties. In the meantime, we may be seriously underesti- 
mating the mass lost in superwinds, and therefore their influ- 
ence within and around the host galaxy. 
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Table 1 

X-ray Components in the Circumnuclear Region of NGC 3079 



Spectral 


n„ a 


kT 


Abundance 


EI b 


Flux c 


Luminosity^ 


region 


(xlO 22 atoms cirT 2 ) 


(keV) 


(relative to solar) 


-3 

cm 


(erg s^'cnf 2 ) 


(erg) 


Nucleus+jet 


1.7±0.2 


0.55±0.05 


0.52 


2.2 ±0.7 x 10 66 


3 x 10~" 


1.1 ±0.1 x 10 42 


Fe line 










9xl0~ 12 


3.3 ±0.5 x 10 41 


Bubble c 


0.3 


0.8 


0.16 


10 66 


8x 10"" 


3 x 10 42 



"Photoelectric absorption using the formulation of Morrison & McCammon (1983). 

b Model component normalizations, where the fitted emission "integral" EI = n e nHrjxV and r\x is the filling factor of volume V. 
Integrated across passband 0. 1-6.5 keV. 
d Distance 17.3 Mpc and isotropic emission. 

e The lack of distinct features in this spectrum prohibits error estimates for this component, and formal \ 2 results in general for 
all components. Where reported, ranges on derived parameters span those values that produced nearly identical "reasonable" fits. 



X-ray/Ha filaments in NGC 3079 
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15"= 1.26 kpc 



Chandra 




FIG. 1 . — (a) Dithered HST/WFPC2 Hq+[N II] A6583 minus I-band image of the superbubble. The "+" spans the current uncertainty in the nucleus position. The 
bubble is formed from vertical towers of filaments anchored to the galaxy disk, numbered 1—4 from south to north. The nuclear wind appears to have "scoured" 
disk gas and dust from within the white dashed circle, (b) 26.6 ksec Chandra/ACIS image at the same scale and orientation, smoothed to 0"8 FWHM. Some regions 
of correspondence A -C are detailed in panel (c); each box is 545 X 360 pc (6"5x4f'3). Ovals on the optical image in A register X-ray features. 




FIG. 2. — Chandra/ACIS spectra extracted from elliptical regions (see text) that are centered on the a) nucleus+jet, and b) superbubble. Solid lines show 
representative fits assuming absorbed, single-temperature models; the nucleus+jet spectrum also uses a single Gaussian to fit the Fe Ka line complex at 6 keV. Fit 
parameters and uncertainties are summarized in Table 1 . 



8 



Cecil, Bland-Hawthorn, & Veilleux 




+55M0'00" 20" 40" +55" , 41'00" 20" 

FIG. 3 . — Large scale X-ray grayscale and radio contours (3.8 cm VLA image from CBVF) are compared. The grayscale is smoothed to 2" FWHM until it reverses 
to show X-ray filaments with 1" smoothing. Contours, suppressed near and immediately west of the nucleus, are drawn at (0.6, 0.55, 0.5, 0.45, 0.4, 0.35, 0.3, 0.25, 
0.2, 0.15, 0.1, 0.05) percent of the peak brightness of 96.2 mJy beam (HPBW 1"6); the lowest two are omitted east of the nucleus to minimize clutter. Also shown 
are the magnetic field vectors. The almost vertical white dotted and dashed lines from CBVF show the modeled positions of the contact discontinuity/ISM shock 
for power-law and exponentially declining gas density above the galaxy disk. They coincide with regions of diffuse X-ray emission and with several high-angle Ha 
filaments shown in gray. 



